Non-damping magnetization oscillations in a single-domain ferromagnet 
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Non-damped oscillations of the magnetization vector of a ferromagnetic system subject to a spin 
polarized current and an external magnetic field are studied theoretically by solving the Landau- 
Lifshitz-Gilbert equation. It is shown that the frequency and amplitude of such oscillations can be 
controlled by means of an applied magnetic field and a spin current. The possibility of injection of 
the oscillating spin current into a non-magnetic system is also discussed. 
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Spin polarized current incident on a magnetic sys- 
tem can exert a torque on its magnetic moment. This 
torque, in turn, can change a magnetic state of the sys- 
tem. One possibility is a switching from a certain mag- 
netic configuration to another one, as predicted theoret- 
icallyl and also observed experimentally in several spin 
valve structures. ^•'^•^ The phenomenon of current-induced 
magnetic switching is a consequence of the spin transfer 
from the conduction electron system to the localized spin 
momentsi^ In certain circumstances, however, the spin 
transfer torque can induce some stable non-damped pre- 
cessional modes. In such states, the energy is pumped 
from the spin current to the localized moments, which 
support the magnetization precession. Such non-damped 
precessions are of high importance from the point of view 
of possible applications in the microwave generation^f^ 

Another important issue in spintronics is the spin injec- 
tion from ferromagnetic to nonmagnetic metals (and/or 
semiconductor) and the spin control over distances com- 
parable to the spin diffusion lengthAiStii Materials 
which might be promising for the applications in spin- 
tronics should have a relatively long spin diffusion length 
(of the order of the system size) and also should allow an 
efficient spin injection across interfaces. Despite several 
technological and fundamental problems, there is some 
progress concerning the efficiency of the spin injection 
and its control by some external parameters. 

From the physical point of view, the phenomena of the 
spin transfer torque and the spin injection are not inde- 
pendent. This is because the current-induced switching 
relies on the spin coherence between two magnetic bodies 
across a nonmagnetic spacing material. In this paper we 
consider precessional modes of a ferromagnetic system, 
driven by a spin polarized current, and the associated 
injection of circularly polarized electrons into a nonmag- 
netic system. This is an extension of our earlier work, 
where we have studied equilibrium and stationary states 
of such a systemii^ To study the magnetic dynamics of 



the system we used the Landau-Lifshitz-Gilbert equation, 
with the spin transfer torque included. We also assumed 
that the torque is an interfacial effect, i.e., the compo- 
nent of the spin current perpendicular to magnetization 
is absorbed at the very interfaceii^ii^ii^ Here we use the 
same model and description to study stable precessional 
modes. 

The time variation of the angular variables 9 and ip, 
which characterize an orientation of the unit vector m 
along the magnetization of the ferromagnetic system, can 
be written in dimensionless variables asi^ii^ 
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with the dimensionless time defined as r = t/[(l -f 
a^)/jHk)], where 7 is the gyromagnetic ratio, Hk is 
the anisotropy field, and a is the damping coefficient. 
Apart from this, h, hs, and hp are the dimensionless ex- 
ternal magnetic field, spin current, and the ea sy-p lane 
anisotropy field, respectively, defined as in Ref. hj. Fi- 
nally, Z{^p) = \ -\- hp cos^ (f. Equations J^l are written 
for the case when the magnetic field and spin current are 
collinear with the easy axis of the ferromagnetic system. 

In the stationary case (r — » 00), the system can 
be transformed into a set of two trigonometric equations 
with respect to the angles and ip^, with the general 
periodic solutions: 



sin2(y9o — ~2hs/hp, cos 6*0 = —h/Z{ip()). 
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The invariance properties for Eqs. ^ allow to select only 
two independent stationary states of Eqs. which can 
be written in the formic 



ifioi + 'P02 = 7r/2, 



coscpo2 = -sgn/is 



Ah^Jhl) 



1/2 



(3) 



2 



sm (^02 



1/2 



sin ipoi = cos (^02 , cos <^oi = sin (^02 , 
cos 6*01,2 -/V2'(<Poi,2)- 

As follows from Eqs. Q, both stationary states and 
their energies do not depend on the damping coefficient 
a, and are determined by ft., hp and hs only. Moreover, 
the latter parameters satisfy the condition < 0.5/ip. 

To investigate the stability of the stationary solutions 
©, we subject them to a small time-dependent pertur- 
bation 

e = 0o + 5ee-"^\ if = ifio + Sip e^"^^ , (4) 
69 — const ^6*0, 6(p = const ^ ipo, 

with uj being generally a complex variable, 

to — UJr + i Aw, (5) 

where ujr is the frequency of homogeneous precession of 
the magnetic moment in an intrinsic effective magnetic 
field, whereas Auj describes damping (or growth) of the 
soft mode fluctuation amplitude and corresponds to the 
natural width of ferromagnetic resonance band. 

Substituting Eq. Q into Eqs. ^ and applying the 
standard linearization procedurei^ with respect to the 
perturbation, one obtains the characteristic equation of 
the second order in u, which - on taking into account 
Eq. ||SJ) - yields the following solutions depending on the 
sign of A = (flu - 022)^ + 4ai2a2i: 



a) A > 
b) A<0: 



= 0, Aw = 0.5 (an + a22 ± Va) , (6) 
±0.5%/^, Aw = 0.5 (an + a22) • (7) 



The matrix elements aij in Eqs. © and Q are defined 
by the stationary solutions Q as 



an = aZ{tfQ) sin 6*0, a2i = Z{lpo) sin 6*0, 
ai2 = hp sin 9q {a cos Oq sin 2ifQ — cos 2(^o) , 
0-22 = hp (cos ^0 sin 2iy9o + Oi cos 2(y9o) ■ 
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As follows from Eqs. © and 0, the oscillatory states 
in the system are possible only for A < 0. The propaga- 
tion of the non-damped oscillations of the magnetization 
components is possible when the conditions A < and 
On = — a22 (Aw = in Q) are simultaneously obeyed, 
which takes place for /is (Aw = 0) = fts^, corresponding 
to the solution of the equation 
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FIG. 1: \uJr\ (surface 1) and 30x |Aw| (surface 2) as a function 
of h and kg for the stationary state (^01, ^01) ■ 



Equation Q is invariant with respect to the simul- 
taneous change of the sign of ft and hsb, h —ft and 
hsb — > —hsb, which corresponds to one of the invari- 
ance properties of Eq. Our analysis has shown that 

Eq. © has one real physical solution. 

The derived formulas to H10|) allow to calculate a 
characteristic dependence of the non-damped oscillations 
of the magnetization vector on the control parameters a, 
ftp, hs and ft (in the linear approximation regarding the 
perturbation). In our calculations we assumed a=0.005 
and ftp=5 (as in Refs. |12m13 ]). reducing in this way the 
number of control parameters to ft and fts . 

Let us consider now the behavior of |wr| and |Aw| with 
ft and fts (Fig. 1). For a better presentation, the values 
of I Aw I were multiplied by a factor of 30 in Fig. 1. As 
one can see, Wr changes slightly with fts and depends 
quadratically on ft, reaching a maximum value at ft = 
0. The surface Aw — f{h,hs), in turn, depends in a 
more complex way on both ft and hg. It is important 
to note that Aw may take either positive or negative 
values, which correspond to an increase or decrease in 
time of the perturbation amplitude. When Aw = 0, the 
system is turned to the neutral mode, when non-damped 
oscillations of a constant amplitude propagate through 
the system. As follows from our calculations, decrease of 
the planar anisotropy hp leads to a further complication 
of the surfaces w^ = f{h,hs) and Aw = f(h,hs), while 
increase of ftp makes them smoother. The increase of the 
damping coefficient a has practically no influence on the 
— fih,hs) surface. 

The non-damped oscillations in th e sy stem are possible 
only when 1 < ft < 1 + ftp (see Ref. [12 )■ This condition 
sets limits on the value of ftsf, as a solution of Eq. I^. 
Figure 2 presents a surface hsb = f{h, hp) corresponding 
to the parameters at which the system is in a state of 
non-damped magnetization oscillations, and thus being 
the boundary between the spin-current stable (fts > hsb) 
and unstable (fts < hsb) states. As follows from Fig. 2, 
hsb depends in a non-linear way on both arguments. 

Consider now numerical solutions of the nonlinear 
equations (1). Figure 3 presents Hausdorff dimensionii 
diagram Z?if(ft, fts) showing main dynamical modes of 
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FIG. 2: Spin current corresponding to the neutral, calculated 
as a function of h and hp. 




Magnetic field h [d.u ] 




FIG. 3: (Color online) Hausdorff dimension diagram 
Dnihjhs) for the system and the main phase portrait types 
for different magnetic fields h and spin currents hs'. 1) h = 
0.4, h, = 0.04; 2) h = 2.825, h, = 0.02; 3) ^ = 12.1, h, = 0; 
4) h = 1.1, hs = -0.02; 5) h = 4.125, h^ = -0.02; 6) 
h = 4.5, hs = -0.02; 7) h = 1.075, hs = -0.03; 8) 
h = 2.925, hs = -0.05; 9) /i = 14.25, hs = -0.03. Dynamic 
modes A ~ F are described in the text. 



the system (phase states) for the different applied fields 
and spin currents. The phase portraits corresponding 
to the most characteristic points of the parameter space 
are shown above and below the bifurcation diagram and 
numbered from 1 to 9. The areas A, B and C correspond 
to the dynamic modes, for which the transition from the 
initial ground state ruz = — 1 to the state = +1 



takes place. In the area A the precession of the magne- 
tization vector takes place mainly along the axis x, and 
the phase trajectory does not cover all the unit sphere 
(phase portrait 1). On the contrary, the phase portraits 

2 and 3, characteristic to the areas B and C, illustrate 
the transition of the phase point to the upper pole via spi- 
ral trajectory, with different initial behavior of the phase 
point running along the two-loop curve (area B) or mov- 
ing from the lower pole along the spiral (area C). Under 
a constant spin current and increasing magnetic field, the 
amplitude of the two-loop curve decreases and the phase 
portrait 2 (area B) turns into that of the phase portrait 

3 (area C). 

When the spin current decreases, the precession of the 
magnetization vector slows down and the phase point be- 
comes unable to reach mz=+l, remaining in the vicinity 
of the two-loop curve and becoming a limit cycle (points 

4 and 5 in Fig. 3, below the upper dashed line and the 
corresponding phase portraits) for the magnetic fields 
1 < h < 1 + hp. Under a further decrease of the spin 
current, the limit cycle of the system turns into a single- 
loop curve (phase 6 and 7), whose form and amplitude 
depend on the control parameters h and hg ■ With a fur- 
ther decrease of the spin current (area D), the limit cycle 
becomes unstable and the magnetization vector instead 
of periodic non-damped oscillations relaxes to a certain 
state with negative and zero rrix (phase portrait 8). 
Upon approaching hum = 1 + hp, the non-stable cycle 
shrinks down to = — 1. It is worth noting that the 
boundary between the areas B and D is well-defined and 
sharp, contrary to the gradual transition between the ar- 
eas A — C. 

When the spin current decreases, the phase portraits 
corresponding to the area C keep the same oscillation 
type, but the magnetization precession becomes signif- 
icantly slower and the resulting spiral trajectory covers 
only a part of the unit sphere (phase portrait 9, area E). 
For the spin currents corresponding to the area F (for 
h < 1 and h > l-\-hp), the phase point is unable to leave 
the ground state niz = —1. Thus, the magnetization 
vector can perform non-damped oscillations, forming the 
phase portraits of a closed cycle for the narrow band of 
the control parameter values (between the dashed curves 
around the boundary between B and D in Fig. 3). Some 
examples of the time dependence of , illustrating non- 
damped oscillations for the phase portraits 5, 6 and 7 
are presented in Fig. 4. The curves in Fig. 4 have been 
plotted starting from a certain time t to eliminate the 
influence of the transition processes taking place in the 
vicinity of r=0. As can be seen in Fig. 4(a), for a given 
hs and increasing h one obtains oscillations of the lower 
amplitude and frequency. Increase of hg (Fig. 4(b), solid 
line) leads to the high-amplitude oscillation with smaller 
frequency. This indicates the possibility of controlling 
the period and amplitude of the non-damped oscillations 
of the magnetization component rriz by the magnetic field 
and spin current. 

To make a qualitative description of the niz (t) behav- 
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FIG. 4: Time dependence of mz{T) for different magnetic 
fields h and spin currents hs'. a) /is = —0.02, solid line h — 
4.125, dashed line h = 4.5; b) hs = -0.03, solid line h = 1.075, 
dashed line - calculations according to 11311 . 



Equation H13|l allows to perform a qualitative descrip- 
tion (in the first approximation) of the possible injection 
of non-damped iriz (r) oscillations from ferromagnetic (F) 
into non-magnetic (TV) material. In the case of an ideal 
injecting contact at x—0, which does not change the value 
and orientation of the spin, one may assume that the 
spin currents to the left and to the right of the contact 
are equalfi^ hf = = /is(0,r). In such a case, the 
functional dependence hs{x, r) can be obtained from the 
continuity equationiS 
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Introducing (|13|l into H14|l . one can show that 
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ior (Fig. 4) we use the expression for 9{t) from Ref. [T^ . 
obtained in the linear approximation with respect to the 
perturbation, 

9{t) = tt + 27rexp(Aa;) sm{uOrT). (11) 

When writing Eq. Hll|l . the initial condition 0{Q) — tt was 
assumed. For the non-damped oscillation mode {Auj = 
0), Eq. (HU yields 

mzir) = — cos [27r sin((jjr''')] . (12) 

As one may expect, the latter equation, obtained in 
the linear approximation, may not give good quantitative 
description of the essentially nonlinear nizir) behavior, 
presented in Fig. 4 as a result of numerical calculations. 
However, one can obtain much better agreement by tak- 
ing the function of the form 

to^(t) = -cos{27rsin^ [^^(t + to)]} , (13) 

where Uro and tq are some approximation parameters. 
For example, the mzir) curve calculated according to 
(HSJ for h, = -0.03, ujro = 1-553 and tq = -0.053 
(Fig. 4(b), dashed line) shows good agreement with 
the corresponding curve obtained by numerical methods 
(Fig. 4(b), solid line). 



hs{x,T) = /is(0,t) -f 7rwrt('?+ -^-) X, (15) 
where f± — cos [2 {tt sin^ LOrQ{T + tq) ± ujra{T + '''o)}] • 

As follows from Eq. (|15|l . the spin current injected from 
the ferromagnetic into the nonmagnetic system will pre- 
serve its non-damped oscillation character, being a su- 
perposition of the second harmonics of two harmonic os- 
cillations. In the framework of the current assumptions 
it will change linearly in amplitude with distance from 
the contact. It is worth noting that the inclusion of re- 
laxation iteml- does not change Eq. (|15|) qualitatively. 
For resistive and other contact typesiiiiSiSS, the expres- 
sion for the spin current will differ from that given by 
Eq. (|15|l . The results of this paper, however, show that 
the non-damped oscillations of mz (t) component can be 
injected from the ferromagnetic to a nonmagnetic system 
due to current continuity at the contact. 



Acknowledgments 

This work is partly supported by FCT Grant 
POCI/FIS/58746/2004 (Portugal), EU RTN2-2001- 
00440 'Spintronics', and Centre of Excellence. V.D. 
thanks the Calouste Gulbenkian Foundation in Portugal 
for support . 



^ J. C. Slonczewski, J. Magn. Magn. Mater. 159, LI (1996); 

195, L261 (1999); L. Berger, Phys. Rev B 54, 9353 (1996). 
^ J. A. Katine, et al, Phys. Rev. Lett. 84, 3149 (2000). 
3 J. GroUier, et al, Appl. Phys. Lett. 78, 3663 (2001). 
* M. AIHaJDarwish, et al, Phys. Rev. Lett., 93, 157203 

(2004). 

^ M. D. Stiles and A. Zangwill, Phys. Rev. B 66, 014407 

(2002); J. Appl. Phys. 91, 6812 (2002). 
^ S. 1. Kiselev, et al, Nature 425, 380 (2003). 

I. N. Krivorotov, et al, Science 307, 228 (2005) 
® J. Xiao, A. Zangwill, and M. D. Stiles, Phys. Rev. B 72, 

14446 (2005). 



^ T. Dietl, Semicond. Sci. Technol. 17, 377 (2002). 
° D. Ferrand, et al, Sol. State Communs. 119, 237 (2001). 
^ E. I. Rashba, Phys. Rev. B 62, R16267 (2000). 
2 P. M. Gorley, et al, cond-mat/0508280 (2005). 
^ J. Z. Sun, Phys. Rev. B 62, 570 (2000). 
" A. Brataas, Yu. V. Nazarov, G. E. W. Bauer, Eur. Phys. 
J. B 22, 99, (2001). 

J. Barnas, et al, Phys. Rev. B 72, 024426 (2005). 
^ N. P. Erugin, I. Z. Shtokalo et al. Ordinary Differential 

Equations (Vyshcha Shkola, Kiev, 1974). 
^ H. Haken, Advanced Synergetics (Springer, Berlin, 1983). 
* J. F. Gregg, et al, J. Phys. D 35, R121 (2002). 



G. Schmidt, et al, Phys. Rev. B 62, R4790 (2000). 

J. D. Albrecht and D. L. Smith, Phys. Rev. B 68, 035340 



